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Abstract. Radiative feedback and magnetic field are understood to have a strong impact on 
the protostellar collapse. We present high resolution numerical calculations of the collapse of a 
1 Mq dense core in solid body rotation, including both radiative transfer and magnetic field. 
Using typical parameters for low-mass cores, we study thoroughly the effect of radiative transfer 
Oh' an d magnetic field on the first core formation and fragmentation. We show that including the 
two aforementioned physical processes does not correspond to the simple picture of adding them 
separately. The interplay between the two is extremely strong, via the magnetic braking and 
' the radiation from the accretion shock. 

a : 



in 



1. Introduction 

The protostellar collapse of low mass dense cores follows a well-defined sequence of 
different stages, down to the formation of a protostar. In the first collapse phase (Larson 
■ 1969), the compressed gas cools efficiently thanks to the coupling between gas and dust. 
At higher densities (p > 10~ 13 g cm -3 ), the radiation is trapped and the first hydro- 
static core (the first Larson core) is formed. At this stage, the grain opacities and the 
radiation transport play a major role. In the recent past few years, a lot of progress in 
the computational star formation field has be done. For instance, a lot of radiation hy- 
drodynamics (RHD) methods have been developed for grid based codes (e.g., Krumholz 
et al. 2007, Kuiper et al. 2010) and for smoothed particles hydrodynamics (SPH) codes 
(Whitehouse & Bate 2006, Stamatellos et al. 2007). Applying these methods to star 
formation, it turns out that a barotropic EOS cannot account for realistic cooling and 
heating of the gas (e.g., Commergon et al. 2010). On larger scales, radiative transfer has 
been found to efficiently reduce the fragmentation thanks to radiative feedback due to 
the accretion and the protostellar evolution (Bate 2009, Offner et al. 2009). Regarding 
magnetic field in the star formation context, a gradually improved expertise has been 
developed for magnetohydrodynamical (MHD) flows (e.g.. Hennebelle Hz Teyssier 2008, 
Machid a et al. 2008|) . All these studies showed that magnetic fields reduce efficiently the 
fragmentation of prestellar cores. Recently, it has been shown that both radiative transfer 
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and ideal MHD arc important and cannot be neglected in the collapse and fragmentation 
of protostellar cores (Price & Bate 2009, Commergon et al. 2010. lTomida et al. 2010[) . 

In this study, we present radiation magnetohydrodynamics (RMHD) calculations of 
prestellar dense core collapse, using the adaptive mesh refinement (AMR) code RAMSES 
(Teyssier 2002). The paper is organized as follows: in section 2, we briefly introduce the 
RMHD solver we designed and present our initial conditions. In section 3, we present our 
results of dense core collapse calculations using the RMHD solver with various numerical 
resolutions. Finally, we draw our conclusion in section 4. 



2. Numerical method and initial conditions 

We use the RAMSES code (Teyssier 2002), in which we have implemented a RHD solver 
using the grey flux-limited diffusion (FLD, e.g. Minerbo 1978) approximation. We use 
the comoving frame, which is valid to leading order in v/c in the static diffusion and 
streaming limits. The solver we designed is coupled to the ideal magneto hydrodynamics 
one developed by Fromang et al. (2006). The coupled RMHD equations read 

d t p + V[/?u] = 

d t pu + V[pu®u-B®B + PI] = -pV$ - AVP r 
d t E T + V [u (E T + P tot ) - B(B.u)] = -pu ■ V3> - P r V : u - AuVP r 

+ V ' (^ V ^) (2.1) 
d t E r + V[uP r ] = -P r V:u + V-(^V£ r ) 

+Kppc(aRT 4 — Er) 
3 ( B - Vx(uxB) = 0, 

where p is the density, u is the velocity vector, B is the magnetic field vector, P to t is the 

"p> TD 

total pressure, sum of the thermal and magnetic pressures, P to t = P H j— , E x is the 

radiative energy, A is the flux limiter (Minerbo 1978), Et is the total fluid energy per 

_ nil — T-l T-l 

unit volume, Et — p(e + ^%^J + E T + 2 , $ is the gravitational potential, P r is the 
radiative pressure, and and up are the Rosseland and Planck mean opacities. This 
system of equations is closed by the perfect gas equation of state P/p = (7 — l)e, where 
7 = 5/3 is the adiabatic exponent. An additional constraint comes from the divergence of 
the magnetic field, which has to vanish everywhere at all times (V.B = 0). Note that the 
radiation diffusion and coupling terms are integrated implicitly in time, since it involves 
very short timescale processes, compared to the hydrodynamical evolution. 

- Initial conditions: We consider a uniform-density sphere of molecular gas, rotating 
about the z-axis with a uniform angular velocity. The prestellar core mass is fixed at 
Mq = 1 Mq and the temperature at 11 K, which corresponds to an isothermal sound 
speed c s o ~ 0.19 km s _1 . To promote fragmentation, we add an m — 2 azimuthal density 
perturbation with an amplitude of 10%. The magnetic field is initially uniform and 
parallel to the rotation axis. The strength of the magnetic field is expressed in terms of 
the mass-to-flux to critical mass-to-flux ratio p, = 20 = (Mo/3>)/(Mo/3>) c - The initial 
ratio of the thermal to gravitational energies is a = 0.37, and the initial ratio of the 
rotational to gravitational energies is f3 = 0.045. 

Calculations were performed using either the rather diffusive Lax Friedrich (LF) Rie- 
mann solver or the more accurate HLLD Riemann solver ( Miy oshi fc Kusano 2005[ ) . Fol- 
lowing up on former studies ([Commergo n et al. 2008[ ), we impose at least 15 cells per 
Jeans length as a grid refinement criterion (parameter Nj). We also use one additional 
parameter, iV exp , which indicates the number of cells refined in each direction around a 
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Figure 1. Temperature and density maps in the yz-plane at time t — l.lhtg for two 
calculations with the FLD and the LF (left) and the HLLD (right) Riemann solvers. 

cell violating the Jeans length criterion. The initial resolution of the grid contains 64 3 



cells. We use the low temperature grey opacities of Semenov et al. (2003) 



3. Results 

- Effect of the solver: Figure [1] shows temperature and density maps in the j/z-planc 
for two calculations with the FLD and the LF and the HLLD Riemann solvers. Each 
calculation has been performed using Nj = 15 and N exp = 4. The FLD-LF case leads 
to spurious fragmentation as it is shown in Commergon et al. (2010) The bubble (dense 
region, p > 10~ 15 g cm -3 ), driven by magnetic pressure due to the magnetic field line 
wrapping, is less extended in the FLD-LF case, and the disc is thus more massive and 
more prone to fragmentation. We identify the accretion shock on the first Larson core 
as a supercritical radiative shock, i.e. all the infalling kinetic energy is radiated away. 
Consequently, the radiative feedback due to the accretion on the first Larson core is 
much larger in the FLD-HLLD case, since the magnetic braking and thus the infall 
velocity are larger thanks to the less diffusive HLLD Riemann solver (Commergon et al. 
2010). 

- Effect of the numerical resolution: Figure shows temperature and density maps 
in the xy-pl&iic for three calculations using the FLD and the LF Riemann solver, with 
various resolutions: (TVj = 15; 7V oxp = 2), (Nj = 15; iV oxp = 4) and (JVj = 20; N cxp = 4). 
We clearly see that increasing the resolution, from left to right, leads to a decrease in the 
number of fragments produced. As resolution increases, the diffusivity of the LF solver 
is reduced, the disk is then less massive, and the magnetic braking more efficient. Using 
the HLLD solver and barotropic calculations, |Commergon et al. (2010) show that the 
correct behavior is the case without fragmentation. 
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Figure 2. Temperature and density maps in the zy-plane at time t — 1.15tff for three 
calculations using the FLD and the LF Riemann solver, with various resolutions. 

4. Conclusion 

We show that taking into account both radiative transfer and magnetic field is not 
a straightforward linear process. We show that the magnetic braking and the magnetic 
bubble extent influence: i) the radiative feedback via the infall velocity and ii) the frag- 
mentation via the disc mass and the rotational velocity. Last but not least, the results 
are extremely sensitive to the numerical resolution and to the numerical diffusivity of the 
code used, which readers and authors should be aware of. 
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